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1. Introduction 

A crucial prediction of Einstein's theory of General Relativity is the existence of 
gravitational waves (GWs), e.g. resulting from the inspiralling and merging process of 
two compact objects. Up till now those waves are purely theoretical predictions with 
lack of direct experimental verification, but their direct detection is under preparation 
by gravitational wave observatories on Earth, e.g., LIGO, VIRGO, GEO, or LISA, a 
future space mission [Tj. 

A first indirect evidence for the existence of GWs was the observation of energy 
loss in the orbital motion of the Hulsc- Taylor binary pulsar PSR B1913+16 being 
in full agreement with the predictions of Einstein's theory using the quadrupolc 
radiation formula. This discovery was awarded the Nobel Prize in 1993. In the 
meantime another strong indirect evidence occurred with the double pulsar system 
PSR J0737-3039A and B |2j [3]. For analysis of the measured GW patterns one has 
to provide very accurate templates following from theory. This can be achieved by 
numerical calculations with the matching of functions to the results or directly using 
analytic tools. In the latter case, waveforms can be obtained only perturbatively due 
to missing analytic solutions of the Einstein field equations for two or more (spinning) 
compact objects (black holes (BHs), neutron stars (NHs)). One of the most successful 
approximation methods is the post-Newtonian one in which the metric keeps close to 
the flat spacetime relying on the assumption that the typical velocity v in a system 
divided by the speed of light c is always small, v/c ~ e <C 1 . The deviations from the 
flat metric can be characterized by the Newtonian potential $; so for a binary system, 
&/c 2 ~ v 2 /c 2 ~ e 2 . In an appropriate limit (as e — > 0), the post-Newtonian (PN) 
approximation yields Newton's equations. 

The merging process of two compact objects is divided into four time scale sectors: 
inspiral, plunge, merger, and ring-down. Each sector delivers characteristic theoretical 
wave patterns, that are hoped to be matched against measured signals in the future. 
The PN approximation provides an excellent analytic handling for the inspiral phase. 
If the PN calculations are very accurate and thus high in order one can make very 
sensible predictions when comparing with measured signals. In this article we focus on 
calculations of the next-to-lcading order (NLO) dynamics of spin-induced quadrupolar 
deformation effects. Surely, the most compact dynamical object is the Hamiltonian, 
generating the equations of motion, so we calculate in section [2] the NLO spin-squared 
one including a constant Cq parameterizing spin-induced quadrupolar deformation 
effects. Cq can be given definite values describing black holes (BHs) or neutron stars 
(NSs). For neutron stars Cq also depends on the model or equation of state (EoS). 
Thus our result, as it seems to be necessary to accurately measure Cq within future 
GW astronomy, could help to find the right EoS. The Hamiltonian in the present 
paper is calculated within the canonical formalism of Arnowitt, Deser, and Misner 
(ADM) [4] . It should be noted that our Hamiltonian is fully reduced in the sense that 
the spin supplementary condition (SSC) is eliminated on the level of the Hamiltonian. 
Further we make a formal counting of the spin as c° and do not distinguish between 
fast and slowly spinning objects (see, e.g., [5] and also Appendix A of [6]). 

By now there are a lot of results regarding spin effects at the conservative orders 
in the PN approximation. The leading order (LO) spin effects are well-known for black 
holes, see, e.g., [2 ® ® [TO] QT] . The LO C Q -dependencc is given in [HUE]. The NLO 
spin effects were only tackled recently. The first derivation of the NLO spin-orbit (SO) 
equations of motion (EoM) is given in [13] which became further developed in |14j . 
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both in harmonic gauge. Later, within the ADM canonical formalism, a Hamiltonian 
presentation was achieved [T5] (see also [IS]). The NLO spin(l)-spin(2) dynamics was 
found in |17[ 116] and confirmed by [181 119] . Higher PN orders linear in spin were 
tackled recently in [501 [UlE]- in particular, Ref. [5D] extended the point-mass ADM 
formalism to spinning objects, valid to any order linear in spin. Even Hamiltonians 
of cubic and higher order in spin were obtained for binary black holes (BBHs) [22| 15]. 
Besides quadrupolar deformations induced by proper rotation (spin) and treated in 
the present paper, tidal deformations induced through the gravitational field of the 
other object were also treated, see, e.g., [2"3" 1 |2"31 |2"5] . 

A nonreduced potential (i.e., with SSC not eliminated on the level of the potential) 
corresponding to the result of the present paper was already calculated in [26], however, 
a term relevant for the center-of-mass motion was missing and only found recently [271 . 
A comparison with the result of the present paper is not trivial, if one wants to avoid 
comparing all (rather long) equations of motion; instead it is more efficient to stay on 
the level of the (relatively short) potential. In section[3]wc sketch how to transform the 
potential from [261 127[ into a reduced Hamiltonian where we will find full agreement 
with our result of the present paper. Considering the special case of black holes (or 
Cq = 1), we already succeeded in calculating the Hamiltonian of the present paper in 
|28[ [5j , providing for the first time both the spin and correct center-of-mass dynamics 
in this case. There we were only able to find agreement with |26) in the spin precession 
equation, see [29] (after identifying a sign typo in [26]). With the correction in [27] a 
full comparison can now be attempted. It will be provided in section [3] 

More work needs to be done for an application of the result of the present paper 
to GW astronomy. In order to obtain the NLO radiation field (for the SO case see 
1301 EH) the stress-energy tensor has to include spin-squared corrections. This stress- 
energy tensor arises from the one with a general quadrupole |32j by a spin-squared 
ansatz for the mass-quadrupole, see [28]. Moreover, the NLO spin contribution should 
be of importance for data analysis. In a recent publication [33j it has been shown that 
for maximal spins (aligned with the total orbital angular momentum), the event rates 
are roughly thirty times larger than of those matter systems with anti-aligned spins to 
orbital angular momentum and eight times as large as for non-spinning binaries. So 
especially considering such sources the event rate will increase with the inclusion of 
spin effects. Further, for the creation of templates, it is useful to find a parametrization 
of the orbits by solving the EoM. It is common to describe the conservative dynamics 
in terms of certain orbital elements. Spin precession and dissipative effects can then 
be formulated as secular EoM of the orbital elements. For explicit solutions including 
spin at LO SO see, e.g., [341135]. 

2. The NLO Spin-Squared Hamiltonian 

We start with giving a short overview concerning the calculation of the Hamiltonian 
in question. This calculation is done within the ADM canonical formalism [J] . We use 
units in which 167rG = c = 1, where G is the Newtonian gravitational constant and c 
the velocity of light. Greek indices will run over 0, 1, 2, 3, Latin over 1, 2, 3. For the 
signature of spacetime we choose +2. We employ the following notations: x = (x 4 ) 
(i = 1, 2, 3) denotes a point in the 3-dimensional Euclidean space R 3 endowed with a 
standard Euclidean metric and a scalar product (denoted by a dot). Letters a and b 
are body labels (usually they are set to 1 or 2), so x„ £ R 3 denotes the position of 
the ath point mass. We also define r a := x — x a , r a := |r a |, n Q := r a /r a ; and for 
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a ^ b, r ab := x a - x 6 , r ab := \r ab \, n ab := r ab /r ab ; \ ■ | stands here for the length 
of a vector. The linear momentum vector of the ath body is denoted by p a = (p a j), 
and m a denotes its mass parameter. The usual flat space spin vector of the ath 
body (in local coordinates) is denoted by S a = (S a u)) in correspondence with our 
paper [28] and its associated antisymmetric tensor by ^(^(j) = UjkS a (k) with the 
total antisymmetric e-symbol defined as £123 = 1. We abbreviate <5(x — x a ) by S a . 
The partial differentiation with respect to x % is denoted by di or by a comma, i.e., 
di(j) = <f) t i; the partial differentiation with respect to x\ we denote by d a i- 

Following the ADM canonical formalism, the independent degress of freedom 
of the gravitational field are described by hf? , the transverse-traceless part of 
hij = cjij — Sij (hJT = 0, hf^j = 0), and by conjugate momenta -k^ tt . The needed 
energy and linear momentum density expressions are given by 



^T^n^nu = H m(attor) , (1) 

-7*?;% = n ? iattcr \ (2) 

where 7 = det(^j), r ) lJ is inverse to 5^, n" is a unit timelike normal to hypersurface 
a; = const, and T^ v is the stress-momentum tensor of the matter system. Hereafter, 
we call its constituents the "particles" , but they may well represent neutron stars or 
black holes. This is substantiated by "general relativity's adherence to the strong 
equivalence principle" : black holes and other compact bodies, to some approximation, 
obey the same laws of motion as test bodies; see, e.g., [36]. Also, the analysis of the 
initial-value solutions for black holes shows that as in electromagnctism, where image 
charges are described by delta functions, black holes in full general relativity can be 
represented by "image masses" with delta functions support |37) . It is convenient to 
choose the following four coordinate conditions 
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tt" = , g tj = + hj/ , = ( 1 + g^J ■ (3) 

The standard ADM Hamiltonian (cf. g]) 

H = j> dSi(g ijt i- 9jj,i), (4) 
then becomes, using the Gauss theorem, 

H = - [ d 3 xA<p. (5) 



The integrand A</> = didi(j) can be expressed in terms of x a , p a , S Q , /i^ T and tt^ tt 
using the constraint equations. By expansions of the field equations in powers of G 
and after adopting suitable regularization procedures of integrals involved (see, e.g., 
Ref. [38] and the Appendix in [39] ). one can determine the Hamiltonian. 

Following the procedure outlined in our previous papers [5] 128] , the Hamiltonian 
and the other generators are constructed as to fulfill the Poincare algebra up to 2PN 
order depending on standard canonical variables 

{xi,Pbj} =&%b, (6) 

{$a \ £>a } = e ijkSa ^ , (7) 

with all other brackets being zero. The coefficient equations resulting from this 
procedure will change due to a modified Hamiltonian and CoM vector entering the 
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crucial relation {Gi,H} — P t see Equation (2.4) in [5]. The modification leading to 
spin quadrupolar deformation effects of a general compact object has to be made in the 
leading order spin-squared Hamiltonian labeled as H 52 in equation (2.8) in [5], which 
now has to include a general spin-quadrupolc constant Cq see [12] ; additionally an 
appropriate Sf-CoM vector has to be found. Both can be accomplished by adopting 
the static source expression for T-L m from our paper [28] equation (4) and incorporating 
the Cq constant reading 

^static = — [h 5l ). ij + 8^ 5m " 7 7 7 -P 1 , q S ^3 S ^l 

7^7™"/ , m §u n Si3kSi) , (8) 

/ ,i 

ci = - \C Q , (9) 

meaning Cq = 1 for BH. Symbolic abbreviations in this formula are taken unaltered 
from the original paper thus denoting the same mathematical objects. This means 
that SVj) being given in an Euclidean basis can be related to a spin tensor Sij in a 

coordinate basis with the help of a triad (dreibein) e^y) by Sy = eitk) e j(l) e klm,S( m ). 
The dreibein as a function of the metric is just = ip 2 Sij because the metric can 
be taken as conformally flat, = ip^^iji m our approximation. The mass-quadrupolc 
tensor of object 1, F£ , is given by 

Il j EE fVV^lfcmSlrU + |s?7 <3 ' , (10) 

2S 2 1 =j ik ~/ jl S lij Siki = const. (11) 

The relation to Q % 1 and & 2 is — m\Q 1 ^ and = mfa 2 , so in leading order the 

related quadrupolc-momcnt tensor is just given by 

Q\ j = afa^ - la?^" (12) 

This static source alone is also enough to determine all the G 2 terms (static, 
free of linear momenta) of the Hamiltonian in question. The LO spin-squared 
Hamiltonian and the 5f-CoM vector are calculated via the formulae H = — J d 3 xA(j> 
and Gi = — J d 3 x x 1 A0, respectively, with a post-Newtonian perturbatively expanded 
A0 and %s2 at s *° r tic according to equations (4.14)- (4.16) in [5]. The results are 

T c Q _ G mrm2 n f (S 1 ■ n 12 ) 2 _ S| 

r 3 
12 



Hsf = ^^C Q ( 3 ^_?" - ^ ) , (13) 



and 



G SI = G— 
1 mi 



(Si- 11x2)8! (S^n^) 2 \ , S 2 , 

^1 2 1 3 (^ 2Xl + ^3X2) + — (^4X1 + ^5X2 , 

r 12 r 12 r 12 



(14) 



with coefficients 

= - \ ~ \c Q • (15) 
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The non-static (momenta based) terms of the Hamiltonian will be determined via the 
same ansatzes for the source terms in W n and H™ as in [5] with the LO quadrupolc 
moment (jT2"j) reading 

2 



w 



E 



6=1 



1 



Pi 



A 



— (Pb ' d)Ql J p bl dj + — &l(p b ■ d) 2 - X s pb ■ (a& x d) Q^didj 



m b 



$b , 



(16) 



Qb hPbkdA + KPbidkdi + A 7 (p 6 • d) 5 h d k 



6=1 L 

+ A 4 a2( Pb . a)a, + ^ (a* x 9), f 1 - iQj'^fliJ - ^6, 



>b ■ 



(17) 



Notice that the static term in H m involves the Cq constant which follows from the 
expansion of <|Sj> being the only modification of Eq. (4.11) in [5]. These sources 
allow the calculation of the NLO spin-squared Hamiltonian (with yet undetermined 
coefficients) leading to the same coefficient equations (4.50) - (4.62) in [S] except for 
the test particle terms 02 and /?4 which are just multiplied by the quadrupolc constant 
Cq . These equations have to be matched to the ones resulting from the requirement of 
fulfilling the Poincare algebra which now include the Cq constant and for that reason 
will slighty differ from the equations (3.8) - (3.22) in [5]. The matching procedure then 
fixes all the coefficients left attributing to the source term coefficients the values 



Ai 



An 



7 V 

4 "2^ 



Ao — 7 + ~Cn, A 



1 

12 



Cq 
6 : 



A (1 



2 

Cq 
4 



1 

24' 



A 4 = 0, 



At 



(18) 
(19) 



which agrees with our results obtained in [5] when setting Cq = 1 (terms cubic in a 
are not of interest here and will be dropped) . In view of section [3] we label from now 
on standard canonical variables with a 'hat' specifying its affiliation to the Newton- 
Wigner (NW) SSC except for the momentum which is chosen to be the same for NW 
and covariant SSC due to the existence of canonical degrees of freedom. The 'hatted' 
variables are then called NW variables in the sense that they are standard canonical, 
meaning 

{xi,p bj }=5 i j 5 ab , (20) 
{S^\S^} = e t]k Si k \ (21) 

Pat = Pai , (22) 

all other brackets being zero. Subtleties arising from that definition of NW variables 
are discussed in detail in our Comment |29] . The resulting NLO spin-squared 
Hamiltonian for general compact binaries reads 



tt ADMcan 
fl NLO S? 



G 



m 2 



21 9 



C Q pftSi-nia)' 



\ + \c Q 
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— ~ ~r C Q (Pi ' A i2)(P2 • ni 2 )(Si • n i2 )' 



3-^C Q )( Pl -p 2 )(S 1 .n 12 ) 2 



21 

T 

9 

r 

3 

3 3 



2 + ? C Q ) (P2 ■ ni 2 )(Si • ni 2 )(Si ■ pi) 



3 + \C Q ) (pi-ni 2 ) (Si -Ai 2 ) (Si • p 2 ) 



(Si- P i)(Si-p 2 ) 



3 3 

o - ~a C Q ) (Pi-ni 2 )(p 2 -ni 2 )Sf 



CV 



-2 + 4^)( pi - p2 ) S i 



^(Si-fii,) 2 --?^ 2 
m 1 m 2 V4 4 



G 2 m 2 



1 „ ?™2 

2 + ~C Q + — 

2 mi 



m 2 



1 + 2Cg Si 



-3--C Q --i(l + 6C Q ) ) (Si-n 12 r 



(23) 



for Cq = 1 being in full agreement with the result for BH presented for the hrst time 
in 125]. 



3. Comparison with the NLO spin(l)spin(l) potential 

In order to transform the Routhian R obtained within the Effective Field Theory 
(EFT) approach [35] [37] to a nonreduced Hamiltonian H, we first have to eliminate 
the acceleration term [27] with the help of the Newtonian equations of motion 
(corresponding to a redefinition of the position variables, see [ID])- This generates 
correction terms to the order G 2 in the NLO spin-squared potential. Next, one 
must replace velocities by canonical momenta p a = J^- to get the Hamiltonian by a 
Legendre transformation, i.e., 

H = vi • pi + v 2 • p 2 - R. (24) 

The canonical momentum p! necessary to cover all NLO spin effects explicitly reads 

(a , 1 A , G?«im 2 
Pi = I 1 + - Vi I rnivi H 2" [6vi - 7v 2 - (m 2 • v 2 )n i2 J 

G 

+ — [m 2 (ni2 x Si) + 2mi(ni2 x S 2 )] , (25) 

r 12 

and similarly for particle 2. The Poisson brackets at this stage are 

{x l a ,Pb 3 } =8}8 ab , (26) 

{SW,S«>} = £l3t Sf), (27) 
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{SM°W} =e i , fc 5(°)W, (28) 
{S (o)( i ) )5 (o)0-) } = _ £i . fcS W ) (29) 

and zero otherwise[jj] Notice that these are not yet the reduced or standard canonical 
brackets as S^ is still an independent degree of freedom and was not eliminated 
using the covariant SSC S^u^ = 0. 

It is well known that one has to proceed to Dirac brackets (DBs) if S^^ is going 
to be eliminated from the Hamiltonian H using a SSC, see, e.g., (42j [41] . However, it 
is possible to find new variables x a , pi,j and Sa for which the Dirac brackets take on 
the standard form, 

{xi,Pbj}r>B = 5)8 ab , (30) 

{S^,S^}vB = e l3k Si k \ (31) 

and zero otherwise. These new variables can only be unique up to canonical 
transformations. This freedom allows us to choose pu = pu, as for the flat space 
A possible transition to x l a and S£ then reads, 



case 



Pi[iSi(j)]( k )Pik 



1 



3 Pl 



4mf 



2Gm 2 
m\f 12 



Pi[iSi(j)]( k )Pik 



3G 
mif 12 
2Gm 2 

9 ^9 



P\\iS\{j)\(k)V\k 



G 



Pi[iSi(j)](k)ni 2 (ni2 ■ p 2 ) 



Pi[i5 , i(j)](0 5 'i(fe)(0"i2 + 



2G 



Pi[iS 1 ( J )](i)S 2 (k)(i)fi h 1 



(32) 



2m\ 



PlkS\(i)(k) I 1 



_pL 



G 



m 2 PikS\(i)(k) 



'■12 



3 rl P2kS\{i){k) 

2 G ' 



Gn k 12 {h 12 •p 2 )5'i( l )(fc) 



min 2 



mir 12 



G 



m 2 S 1{k ) {l )S 1{i )(i)h{ 2 i G ,n 12 5 1 ( i )( i )S , 2(fc)(i) 



TOT 



TOif 12 



(33) 

n ' 12 

where the antisymmctrization of indices pertaining to a tensor Ay is defined as 
A 



1/2 ( Aij — Aji). The rather complicated form of these variable transformations 



reflects the complicated structure of the DBs for self-interacting spinning objects in the 
covariant SSC. We will elaborate on its specific calculation in another paper. Notice 
that these results are applicable to all NLO spin effects (for the spin-orbit contributions 
the corrected form of the tetrad from [25] has to be inserted into the SSC). To best 
of our knowledge this is the first time that DBs are applied to gravitationally self- 
interacting spinning objects. In [2T] test-spinning objects are considered and [5T] 
covers the flat space case only. 

The Routhian from [26] now leads us to the reduced NLO spin-squared 
Hamiltonian in the form 



rrEFTcan 



G 



m 2 



21 



c Q \pi(s 1 -h 1 2y 



II Notice that in 11811 261 a different sign convention was used for the Poisson brackets of the spin. As 
usual, we are not showing the canonical conjugate of the spin here; see [41] . 
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2 -i+ 9 -C Q )( Pl .n 12 fSl 
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4 " 

2q2 



Cq pfSf 



2^Q)(Si-Pir 



1/15 

~ 2 ~ l - ^ (Pi ' ft i2)(P2 • ni 2 )(Si • n 12 ) 5 



-3+-C q ) ( Pl -nia)(Si •n 12 )(S 1 • p 2 ) 



2--Cg ) (Si-pi)(Si-p 2 ) 

-2 + 



4 ^ 



(pi • n i2 )(p 2 • ni 2 )Si 



TOlTO 2 \4 4 



G 2 m 2 



2 mi \2 



m 2 (\ 



+ I (-+6C, 



mi \2 



Q 



3Cr 



(Si-n 



Si 



12 J 



(34) 



This Hamiltonian and the one calculated in the last section should differ only up 
to a canonical transformation. It should thus be possible to generate the difference 



AH 



ff^LOSf 1 — ^NLOSf" 1 a canon i ca l transformation of the form 



A-HnlosI = i H N, 9^032} , (35) 

with Hn being the Newtonian Hamiltonian of a two-body system and g being an 
appropriate generator. It turns out that with the generator 



can 
i>NLOS? 



G m 2 



+ C Q (pi-n ia )S? + -(S 1 -n 12 )(Si-pi) 



(36) 



equation (|35[) can be fulfilled and so agreement is achieved. This means that the 
Hamiltonian H^^f n calculated with the aid of the ADM method, in terms of 
invariant physical quantities, agrees with the Routhian from above, hence there is 
great confidence that -ff^LOS? 11 correctly describes a binary consisting of BHs and/or 
NSs or other kinds of compact objects in post-Newtonian Einsteinian theory. The 
new Hamiltonian may find immediate application in the problem of motion of orbiting 
binaries as investigated and solved in e.g., [3SJH3]. 
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